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Assessment of a 3-D Boundary Layer Analysis to Predict Heat Transfer and 

Flow Field in a Turbine Passage 


1.0 SUMMARY 


An assessment has been made of the applicability of a three dimensional 
boundary layer analysis to the calculation of heat transfer, total pressure 
losses, and streamline flow patterns on the surfaces of both stationary and 
rotating turbine passages. In support of this effort, an analysis has been devel- 
oped to calculate a general nonorthogonal surface coordinate system for arbitrary 
three dimensional surfaces and also to calculate the boundary layer edge condi- 
tions for compressible flow using the surface Euler equations and experimental 
pressure distributions. Using available experimental data to calibrate the 
method, calculations are presented for the pressure, endwall, and suction sur- 
faces of a stationary cascade and for the pressure surface of a rotating turbine 
blade. The results strongly indicate that the three dimensional boundary layer 
analysis can give good predictions of the flow field, loss, and heat transfer on 
the pressure, suction, and endwall surface of a gas turbine passage. 
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2.0 INTRODUCTION 


The prediction of the complete flow field in a turbine passage is an 
extremely difficult task due to the complex three-dimensional flow pattern which 
contains separation and attachment lines, a saddle point, and a horseshoe vortex 
(Fig. 1). Whereas, in principal such a problem can be solved using full Navier- 
Stokes equations, in reality methods based on a Navier-Stokes solution procedure 
encounter difficulty in accurately predicting surface quantities, such as heat 
transfer, due to grid limitations imposed by the speed and size of the existing 
computers. On the other hand the overall problem is strongly three dimensional 
and too complex to be analyzed by the current design methods based on inviscid 
and/or viscous strip theories. Thus there is a strong need for local enhancing 
of the current prediction techniques through inclusion of 3-D viscous effects. A 
potentially simple and cost effective way to achieve this goal is to use a 
prediction method based on three dimensional boundary layer (3-DBL) theory. The 
major objective of this study is to assess the applicability of such a 3-DBL 
approach for the prediction of heat loads, boundary layer growth, pressure 
losses, and streamline skewing in critical areas of a turbine passage. For this 
purpose, the three dimensional boundary layer analysis developed by Vatsa (Ref. 1 
and 2) has been selected to evaluate this approach as a means for calculating the 
local properties of the flow field. 

In this approach zonal concepts are utilized to delineate regions of appli- 
cation of 3-DBL theory - these being the endwall surface, suction surface, and 
pressure surface of a turbine blade as shown by the shaded regions of Fig. 1. The 
zonal concept employed in this study implies that there exists a thin region near 
the surface dominated by wall pressure forces, friction forces, and Coriolis 
forces so that boundary layer theory is valid provided that the proper inflow 
conditions and boundary layer edge conditions are specified. Although the pres- 
sure surface of a stationary blade (cascade) shows only weak three dimensional 
effects, the suction surface shows strong effects due the nearby passage vortex 
which sweeps the flow from the endwall. Likewise the pressure surface of a rotat- 
ing turbine blade shows strong three dimensional effects due to the interaction 
of the strong radial pressure gradient and the Coriolis force. These strong three 
dimensional effects should provide a rigorous test of the zonal application of 
3-D boundary layer theory to the turbine. 

This zonal approach requires three separate analyses: 1) an analysis to 
construct a general non-orthogonal surface coordinate system for twisted turbine 
blades, 2) an analysis to calculate the boundary layer edge conditions from a 
known static pressure distribution, and 3) a 3-D boundary layer analysis which 
predicts the boundary layer growth with prescribed inflow conditions. A review 
of the background literature on these three problems is given below. 

A coordinate system must have certain general properties if it is to be 
useful for calculating three dimensional boundary layer solutions. Since the 
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boundary layers lie on the turbine blade surface, a useful coordinate system 
would be one formed by the intersection of three sets one parameter surfaces of 
which two sets of surfaces intersect the wall boundary and the third set of 
surfaces move off the wall boundary in a one parameter set. In such a coordinate 
system, one surface of the third set is the boundary surface which is described 
by only two parameters (coordinates). Thus as an example Howarth (Ref. 3) 
derived the three dimensional boundary layer equations in a general orthogonal 
coordinate system which fits this requirement. However triply orthogonal coordi- 
nate systems are difficult to construct for arbitrary surfaces such as a twisted 
turbine blade. Squire (Ref. 4) derived a more general set of of boundary layer 
equations in a restricted nonorthogonal coordinate system. However this set of 
equations has certain coordinate curvature restrictions which make it difficult 
to apply in practical cases. If, however, one makes the assumption that the 
boundary layers are very thin compared to the radius of curvature of the surface, 
then the problem is greatly simplified. In this situation the third set of 
surfaces is approximated by the boundary surface and is called a surface coordi- 
nate system. In this surface coordinate system, two coordinates lie on the 
surface and the third is normal to the surface. In addition it should be noted 
that on the boundary surface, it is difficult to construct two orthogonal 
families of curves to describe the surface. Thus it is desirable to have a non- 
orthogonal coordinate system on the surface. The problem then reduces to the 
mapping of the three dimensional boundary surface to a plane surface to be 
described by two one parameter families of curves (coordinates). This mapping 
function has been developed by Gordon and Thiel (Ref. 5) and the construction of 
a general nonorthogonal surface coordinate system for arbitrary surfaces is 
described in this report. Since turbine blades are rotating, another requirement 
of the coordinate system is that it should be a rotating coordinate sysem so that 
the Coriolis, or apparent forces, appear explicity in the boundary layer equa- 
tions. Mager (Ref. 6) has derived the boundary layer equations in a general 
orthogonal rotating coordinate system. However, as stated earlier, a nonortho- 
gonal surface coordinate system is more useful. Vatsa (Ref. 1 and 2) has derived 
a set of boundary layer equations in a nonorthogonal rotating surface coordinate 
system which meets all these requirements and therefore this analysis has been 
used to make the assessment presented in this study. 

The solution of the boundary layer equations requires specification of the 
boundary layer edge conditions. These conditions are the two components of the 
edge velocity, the edge total enthalpy (rothalpy), and the thermodynamic vari- 
ables of state. These edge conditions can be obtained directly from experimental 
data or they can be obtained from solutions of the Euler equations. Measurements 
of the vector velocity components in a three dimensional flow field are extremely 
difficult and costly to obtain. Difficulties are also encountered in the use of 
the Euler equations for the boundary layer edge conditions since these solutions 
do not produce secondary flows which are generated by viscous shear forces. An 
alternative approach is to obtain the edge conditions by solving the Euler 
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equations evaluated at the surface (hereafter referred to as surface Euler equa- 
tions) using a known (experimental) pressure distribution. This approach fits 
well within the scope of the present program since the overall objective is the 
assessment of 3DBL analysis for turbine flows. Since static pressure distri- 
butions over surface are relatively easy to obtain, this method is more straight- 
forward and avoids the problems mentioned above. This method was outlined by 
Cebeci (Ref. 7) for application to aircraft wings. Gleyzes and Cousteix (Ref. 8) 
developed a similar method for incompressible flow over fusiform bodies. This 
report extends this method to the prediction of the boundary layer edge veloci- 
ties and thermodynamic quantities for compressible flow over general three dimen- 
sional surfaces. 

A detailed review of the development of three dimensional boundary layer 
theory is given by Vatsa (Refs. 1 and 2). In this report only a brief outline of 
the methods used by Vatsa shall be given. It has long been recognized that 
turbulent boundary layer growth is governed by two length scales that have 
different properties. Near the wall the turbulence is affected by the presence of 
the wall and the inner length scale reflects this property of the turbulence. 

Thus flow is described by the well known law of the wall. Far from the wall, the 
turbulence is wakelike in behavior and an outer layer length scale describes the 
turbulence properties. For two dimensional laminar boundary layers, the Levy-Lees 
transformation, such as that used by Blottner (Ref. 9), attempts to capture the 
growth of the boundary layer and thereby significantly simplifying the analysis. 
For turbulent boundary layers, Werle and Verdon (Ref. 10) have generalized this 
concept by replacing the molecular edge viscosity with an effective turbulent 
viscosity. Vatsa (Ref. 1 and 2) has generalized these concepts to three dimen- 
sional turbulent boundary layers and has sucessfully obtained solutions to a 
number of problems. 

In this report, an assessment of the applicability of a three dimensional 
boundary layer analysis to the calculation of heat transfer, total pressure loss, 
and streamline skewing on turbine blades is made using the 3-D boundary layer 
analysis of Vatsa (Refs. 1 and 2). In support of this assessment an analysis has 
been developed to construct a general nonorthogonal surface coordinate system for 
arbitrary three dimensional surfaces and an analysis has also been developed to 
calculate the boundary layer edge conditions using the surface Euler equations 
and experimental surface static pressure distributions. Both of these analyses 
are explained in detail in Section 3 - Analysis, along with a review of the 3DBL 
analysis developed by Vatsa. In Section 4 - Results and Discussion, available 
experimental data is used to calibrate the method with calculations presented for 
the pressure, endwall, and suction surface of a gas turbine cascade described by 
Graziani et. al. (Ref. 11) and for the the pressure surface of a rotating 
turbine blade described by Dring and Joslyn (Refs. 12 and 13). 
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3.0 ANALYSIS 


3.1 Surface Coordinate System 


The three dimensional boundary layer equations are written in a surface 
coordinate system (x p x 2 , x 3 ) in which x^ and x 2 lie on the surface and x j is 
orthogonal to (xj, x 2 ) and hence normal to the surface. The coordinate Xj is 
generally in the streamwise direction and x 2 is generally in the crossflow 
direction. If the coordinates of the surface are written in Cartesian 
coordinates (y^ y 2 , y 3 )» and the transformation 


Yi - y £ (*j> 


(3.1.1) 


is known, where the Jacobian 


J - |*1| ,0 

a Xj 


(3.1.2) 


then the components of the covariant metric tensor are given by Warsi (Ref. 14) 


g i i 


_ 9y k 9y k 


^ 3x^ 3xj 


(3.1.3) 


For the surface coordinates defined above, Eq. (3.1.3) reduces to 


811 


/jM 2 , /M 2 1 /iaf 
l 3x^ J l 3x^ J l 3x^ J 


(3.1.4) 


g 12 


9y l 9y l + 9y 2 9y 2 + 9y 3 9y 3 


3x^ 3 x 2 


3x^ 3 x 2 


3x^ 3 x 2 


(3.1.5) 


g l 3 


= 0 


(3.1.6) 
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§21 ®12 (3.1.7) 

-■ ■ • (§)’ • ( 5 )' 


823 = 

0 

(3.1.9) 

831 = 

0 

(3.1.10) 

832 = 

0 

(3.1.11) 

833 = 

1 

(3.1.12) 


The determant of the metric tensor is given by 

|g I = 811822 “ ® 12 ^ (3.1.13) 


and the metric scale coefficients are 


hj - /g^ (3.1.14) 

h 2 = ^822 (3.1.15) 


h 



(3.1.16) 


Arc lengths along the coordinates are then determined by the relations 


ds^ = h^dxj 


ds 2 h^dx^ 


ds^ = dx-j 


(3.1.17) 

(3.1.18) 

(3.1.19) 


6 



R85-956834 


It is noted that the angle between x^ and X 2 in nonorthogonal coordinates is 
given by 


cos 9 = g 12 /(h 1 h 2 ) 


(3.1.20) 


The surface coordinates (xp x 2 , x 3 ) are the computational coordinates. The 
Cartesian coordinates (y^ y 2 , y 3 ) , which are used as a basis for the 
transformation, are the physical coordinates. The transformation, Eq. (3.1.1), 
defines a unique point in physical space with its corresponding point in 
computational space. In order to insure uniqueness, the Jacobian Eq. (3.1.2), 
must never pass through zero anywhere in the computational domain. This occurs 
at the leading edge and trailing edge of a turbine blade. Therefore the 
computational domain must extend from a point just downstream from the leading 
edge of the blade to a point just upstream of the trailing edge of the blade. In 
general (xp x 2 , x 3 ) do not represent physical distances along the surface of the 
blade. The physical distances along the coordinates, which lie on the surface 
of the blade, are obtained from the metric scale coefficients using Eqs . (3.1.7 
through 3.1.19). A special case of surface coordinates is orthogonal coordinates 
on a plain surface. For this special case we have 


§11 = 1 

(3.1.20) 

g 12 = 0 

(3.1.21) 

g 22 = 1 

(3.1.22) 


and, as can be seen, the physical and computational coordinates are identical. 
Although it is possible to construct other simple coordinate systems by analytic 
means, these coordinates are not useful for turbine blades. Since turbine blades 
are twisted surfaces and their coordinates are not known except in numerical 
form, special analysis is required to construct a general surface coordinate 
system. This analysis is described in Section 3.10. 


R85-956834 


3.2 Boundary Layer Equations 

An ideal set of equations for this problem consists of the three dimensional 
boundary layer equations in a general nonorthogonal rotating coordinate system so 
the Coroilis forces appear explicitly in equations. In addition to the general- 
ized geometry required for realistic turbine blades, it is noted that boundary 
layers on turbine blades are laminar, transitional, and/or turbulent. Therefore 
the boundary layer equations should include the Reynold's stress components so 
that the appropriate turbulence models can be applied. Vatsa (Refs. 1 and 2) has 
derived a set of three dimensional turbulent boundary layer equations that meet 
these requirements. This derivation is long and involved and therefore will only 
be outlined in this section. The first step is to transform the Navier-Stokes 
equations from a stationary Cartesian coordinate system to a moving coordinate 
system using the Galilian transformation. 


y = y + v b t (3.2.1) 


■* * * 

u=u+v B (3.2.2) 


v B =J)xr (3.2.3) 


The second step is to derive the Reynolds stress terms and include them with the 
molecular viscous stress terms by taking an ensemble average over all possible 
instantaneous flow conditions, 


1 

f = 1 im — E f 

N -*■ » N 


(3.2.4) 


where f represents any combination of dependent variables. All products of the 
dependent variables are then separated into its average and fluctuating compon- 
ents (correlations). 
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fg = f g + f g 


(3.2.5) 


The fluctuating components are then included with the molecular stress terms. 

For compressible flow, the number of terms is greatly reduced by neglecting terms 
with triple correlations and correlations with density fluctuations. The Navier- 
Stokes equations are then transformed from the Cartesian coordinate system (y x , 
y 2 j Y3> to a general nonorthogonal coordinate system (x A , x 2 , x 3 ) and then 
reduced to a surface coordinate system using the simplified covarient metric 
tensor components defined by Eqs. (3.1.4) through (3.1.12). Finally the 
equations are simplified using the boundary layer assumptions. Using this 
procedure Vatsa (Refs. 1 and 2) has derived the following set of equations using 
the tensor relations given by Warsi (Ref. 6). 


Continuity Equation 
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- 2 


h l h 2 


<i> 2 u 2 


g 12 


w 3 u l 


2 h 1^2" 

- arr — — — ■ 



2 

a) r 


hjg 12 3r 

g 9 x 2 


O 

h x h 2 3p h i g 1 2 3 P 

Pg 9x 1 Pg 3 x 2 


1 

P 






(3.2.7) 


X 2 Momentum Equation 
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3r 
g 


2 "2 812 Vl 2 “1 “2 

CD r — = — - tD z r 1 z 


hi h 0 3r 


g 3 Xr 


h 2 g 12 3p h l h 2 3p 

Pg 3xi Pg 3x 9 


1 3 


( 3u 2 

l y — - - p u,'u,' 
p 3x 3 J 3x 3 


Energy Equation 


u j 3hrji U2 3 hfj* 

h^ 3x^ h 2 3 x 2 


3h 


T _ 


3xo 


1 3 


P 3x 3 ( 


i 1 *1 


9x- 


P u' 3 h r ' + y 


3x 0 


(^)i 


U 1 


3 00 u 9 3 

+ (co 2 r 2 ) + — =— 

2h l 3 X j 2h 2 3X 2 


(cd 2 t 2 ) + — ( (D 2 r 2 ) 

2 3xo 


In addition to the equations of motion, we have additional relations 
given below. 

Equation of State 


p = P t 


Stress/Strain Heat Flux Relations 


3ui 


T 13 = M T— - _ P Ui'ut' = (y + £,) 


3x 


1 3 


3uj 

3xo 


3u^ 


T 23 = u “ p u 2'u^' = (y + e 0 ) 


3X' 


2 3 


3u 2 

3Xo 


(3.2.8) 


(3.2.9) 


which are 


(3.2 .10) 


(3.2.11) 


(3.2.12) 
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= X 


at 


- p u 3 'h t ' 




y 

Pr 


a 

ax 3 



(3.2.13) 


where and e 2 are the nonisotropic components of the eddy viscosity, qj is 
the eddy conductivity, and the turbulent Prandtl number is defined by, 


P rt C P e l /e H 


(3.2.14) 


Total Velocity 




+ 2 


u l u 2 


£12 

h l h 2 


(3.2.15) 


where the third term is due to the nonorthogonality of the coordinates. 
Total Enthalpy 


h 


T 


V + 



Rothalpy 


i = hj - Vg^/2 


(3.1.16) 


(3.1.17) 


where the rotor speed is given by 


V B = ru) 


(3.2.18) 


For stationary coordinates, oo = 0 and the rothalpy is identical to the total 
enthalpy. On a rotating blade, rothalpy is conserved along a streamline. On a 
stationary blade, total enthalpy is conserved along a streamline. Except for 
very special cases, it cannot be assumed that the rothalpy is constant over the 
surface of a rotating blade. 
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Sutherland's Viscosity Law 


V 


w ref 


(l \ 3/2 t Mf * 198.6 

l^ref/ ‘ * 198 ' 6 


(3.2.19) 


A number of comments should be made about these equations. First it is 
noted that these equations are valid only for very thin boundary layers relative 
to the radius of curvature of the surface. For thick boundary layers on highly 
curved surfaces a more general coordinate system is required to properly account 
for the divergence of the streamlines at the edge of the boundary layer. 

Secondly, it is noted that certain stress terms related to surface curvature are 
neglected. However, as Bradshaw (Ref. 15) has pointed out, the effects of wall 
curvature on the generation and decay of turbulence are far larger than the 
neglected terms. Therefore curvature effects are best treated in the turbulence 
model. Likewise any effcts of Coriolis forces on the generation and decay of 
turbulence can also be treated in the turbulence model. Thirdly, it is noted 
that only one component of the Coriolis force appears in the boundary layer 
equations. This component is normal to the wall which acts to turn the flow in 
the Xj or x 2 direction. For laminar flow, the equations given above form a 
complete set of equations. For turbulent flow, a turbulence model is required to 
determine the effective viscosity and conductivity of the flow. The turbulence 
models used in this study are treated in Section 3.8. 


3.3 Normalized Equations 

The equations in Section 3.2 are normalized with respect to a reference 
length (jt), a reference velocity (u^), and a reference density (p,,,). The 
reference temperature is defined in terms of the reference velocity and the gas 
constant^. In addition, since the boundary layers are very thin, the coordinate 
Xg and the velocity u^ are scaled by the square root of the Reynolds number. 

Thus all variables are normalized as follows: 


Xj = Xj/£ 

X 2 = x 2 /fc 
X 3 = x 3 /£/Re 


U 

U 

U 


1 

2 

3 


Uj/Uc 

u 2 / u «o 


Ua/u /Re 

3 m 


(3.3.1) 


(3.3.2) 
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2 

P = p/( PeoU^ ) 

p = p/ p„ 

t ref = u « 2/ ^ 

T = ^ref 

% = h t /u M 2 

i = i/ Uoo 2 


ft = u)£/u 

R = r / 1 

y ref = y (t ref ) 

Re = Poo U oo^/y re f 

y - y/u ref 


H 1 = h x /l 
H 2 = h 2 /Z 

G 12 = ® 12/ * 2 
|G | = |g |/* 4 


(3.3.3) 


(3.3.4) 


(3.3.5) 


(3.3.6) 


The resulting normalized equations are the same as those given in Seciton 3.2 
with lower case letters replaced by upper case letters. The equation of state 
and the total enthalpy change as follows; 


P = pT 


(3.3.7) 


H„, 


■y 


t + 



(3.3.8) 
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3.4 Generalized Levy - Lees Transformation 

Boundary layers driven by pressure gradients undergo significant growth or 
contraction which is difficult to estimate apriori. In addition, turbulent boun- 
dary layers have two length scales; an inner length scale near the wall 
reflecting properties of the law of the wall, and an outer length scale 
reflecting wake like turbulence behavior. These properties of the turbulent 
boundary layer make selection of the grid distribution extremely difficult if the 
flow were solved in physical variables. For laminar two dimensional boundary 
layers, the Levy-Lees transformation, as given by Blottner (Ref. 9), effectively 
captures the boundary layer growth thereby significantly simplifying the 
analysis. For turbulent boundary layers, Werle and Verdon (Ref. 10) have 
generalized this concept by replacing the laminar edge viscosity coefficient with 
an effective turbulent viscosity coefficient resulting in a turbulent version of 
the Levy-Lees transformation. For three dimensional boundary layers, Blottner 
(Ref. 16) has reviewed some of the transformations currently used. While these 
transformations work reasonably well for laminar or turbulent flows, they are not 
entirely satisfactory over the complete range of laminar, turbulent, and 
transitional flows encountered on gas turbine blades. Vatsa (Ref. 17) has 
developed a more suitable transformation by extension of the Levy-Lees variables 
to the three dimensional boundary layer equations. This transformation reduces 
to those of Blottner (Ref. 9) and Werle and Verdon (Ref. 10) for two dimensional 
boundary layers. The use of these Levy-Lees variables also has the advantage, 
which will be exploited in Section 3.7, of allowing the calculation of a family 
of similarity solutions to be used as inflow conditions for starting the calcula- 
tion. 


The generalized Levy-Lees transformation uses independent variables defined 


by. 


* 1 = 1 % dX l 


(3.4.1) 


h ~ X 2 


(3.4.2) 


^3 = 


U le H 2 


VTjj- l » dX 3 


(3.4.3) 
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where 


q " p e U le H l H 2 2p e e ref 

(3 

o 

II 

X 

(-• 

o 
^ — ' 

(3 

= 1 + Te/y 

(3 

e re f = 1 + T (e/y) e 

(3 

L = py 

(3 

P e b e 


The subscript (e) is used to denote the boundary layer edge. 
New dependent variables are defined by 


F = Uj/U^ 

(3 

G = U 2 /U ref 

(3 

0 = T/T e 

(3 

H = H T /H Te 

(3 


where U re f may be either or U 2e . In these new variables a transformed 
velocity is defined by integrating the continuity equation (3.2.6), 


— n, n 0 

V = /25 1 q P W U 3W 



H 1 H 2 





H 1 H 2 


a 





J ref 



(3 


4.4) 

4.5) 

4.6) 

4.7) 

4.8) 


4.9) 

4.10) 

4.11) 

4.12) 

normal 


4.13) 
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where is the wall injection velocity. The boundary layer equatioi 
generalized Levy-Lees variables as derived by Vatsa (Ref. 1) are then 
follows : 

Xj Momentum Equation 


H 3 


'ref 


3F 

9£c 


- V 


3F 

3£ 3 


- A 1 F 


3F 3F 

- Aj G — 


H 


1 


H- 


- A 12 f2 " A l 3 gf " A 7 G 2 + AgF + AgG 


+ A 14 ® + A 15 - ® 


where the coefficients are given by 


A 1 = 2 ^1 


A 2 = 2 


Ci Hj U 


ref 


q H 2 U 1( 


q„ 1 3U, 

A 3 = 2 ?1 — — — 

1 q u i e Hi 


A 4 = 


U ref 1 9U le 

q H 2 U le U le 3C 2 


K 1 = 


’12 3H, 3H, q 

q + - — 

2 q ° 


i i 


1 


3 G 


12 


Hi Hj H, 3C 


«i g 12 
A 5 = — H 1 K 1 


s in the 
written as 


(3.4.14) 

(3.4.15) 

(3.4.16) 

(3.4.17) 

( 3 . 4 . 18 ) 

(3.4.19) 

(3.4.20) 
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K 


3 



- H 2<o 


3H 2 


G 12 3H 2 

IT ~^2 


A 


7 




K 


3 


A 8 ~ 




H, 


H, 


'l G l U 


U ref ^3 


4 Cj H 1 G 12 % 

<J S\G\-~ U le 


Ai n ~ 


2 ?1 


H, 2 H 2 3R , 

— r— q o R — 

G| 


(3.4.21) 


(3.4.22) 


(3.4.23) 


(3.4.24) 


(3.4.25) 


(3.4.26) 


(3.4.27) 
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2 q 

A, , = i 

Hi 2 S 12 

3R , 

P 

(3.4.28) 


AH 

q 

M u le 

h 2 



A 12 = A 3 + 

A 5 

(3.4.29) 



A 13 = A 4 + 

A 6 

(3.4.30) 

A 14 " A 3 + 

( a 4 + A 6 “ A 8 ^ 

G 0 + A 7 G 2 
e / e 

+ a 5 ~ A 9 - Ajq - a x1 

(3.4.31) 



A 15 = A 10 + 

A 11 

(3.4.32) 

Momentum Equation 




3 

T "2 9G 

3G 

- y 

3G 3G 

- A, F - A, G 


"9% 

e ref 9? 3 

9 ?3 

1 Hi Z 95 2 



— A 27 G “* ^28 A 23 G 

- a 20 f 2 - a 24 f 



+ A 30 0 + A 29 = 0 (3.4.33) 


where 


A i 6 “ 2 ?i -r 


3U 


2 e 


q U ref 9 ^1 


q 1 

A 17 = 2 ?1 -5 


3U 


ref 


* U ref 9 ?1 


(3.4.34) 


(3.4.35) 
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.2 £l Hj U 2 e 1 ^2e 

q H 2 U ref U le *2 

- 2 ^1 H J_ 1__ l^ref 

q H 2 U le 9 ?2 


9G 17 9H, G 17 9Hi 

K 4 = q o —if- - H, 1 - 12 - 1 


H 


1 


1 H 2 h i ° Hi 


v 20 


2 Ci Hj H 2 U le 


q G | U 


ref 


k 5 = h 1 h 2 


1 + 



9Hj 


- 2 G 


12 


9H] 
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. 2 ^1 H 1 


q g 


g 12 9H 2 1 9G 12 3H 2 


6 u 2 r±r U ^r- % 


2 ^2 n 2 


Hi 


. _ 2 ?1 H 1 G 12 U ref „ 

A 22 - — — -7-T Kg 


q g 


u 


le 


(3.4.36) 


(3.4.37) 


(3.4.38) 


(3.4.39) 


(3.4.40) 


(3.4.41) 


(3.4.42) 


(3.4.43) 
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^23 


^ ^1 ^ 1*^12 % 

q /|G| ^le 


(3.4.44) 


v 24 


4?i ^3 

q 4TT u re f 


(3.4.45) 


v 25 


_ 2 h H 1 G 12 'lo 


CT 3R 

R 


U ref U le H] 


(3.4.46) 
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3R 


l G l U le U ref 9 ?2 


(3.4.47) 


A 27 A 19 + A 22 


(3.4.48) 


A 28 = A 17 + A 21 


(3.4.49) 


A 29 A 26 A 25 


(3.4.50) 


a 30 “ a 16 + a 18 + A 20 + A 24 ~ A 29 + A 21 G e + A 23 G e + A 22 G e (3-4-51) 


Energy Equation 


9 % 




3H 


Pr e. 


ref 9? 3 


- V 


3H 


- AjT 


3H 


- A 2 G 


3H 


Hi H 2 


- A^ FH - A^2 GH + A^F + A^ G 
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9F 

+ 

9 

Aoq LG 

35 3 

9G 

ac 3 

9^3 

9G 

9F ) 


F 

+ LG > 

= 0 

9^3 

3C 3 ) 



% 1 

9H e 

q 

H e 

3?1 

H i 

U ref 

3H e 

[ 2 

U le 

3?2 


3 R 

9R e 

q 

H e 

95 1 

H i 

U ref 

R i 

h 2 

U le 

H e ; 


«i. 



H e 



U ref 


» 

H e 



U ref 

G 12 


U, _ H 


(3.4.56) 


(3.4.57) 


(3.4.58) 


(3.4.59) 
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A 39 “ A 36 A 38 
A 40 = A 37 A 38 


Continuity Equation 


3V 3F 


3G 


+ A 1 — + A 2 TT + A 45 G + A 46 F " 0 


9 % * 9 ?! ^ 9? 2 


where 


A 41 ' 2 


H 1 9 


q H 2 3 ?2 




K ? -M 


H 2 Hi 


42 


_ A i 3 K ? 


(3.4.60) 

(3.4.61) 

(3.4.62) 

(3.4.63) 

(3.4.64) 

(3.4.65) 

(3.4.66) 


K £ 


= ^|G| 


Hr 


*43 


A 2 3 K g 


K, 


8 9? 2 


*44 


(3.4.67) 


(3.4.68) 


(3.4.69) 
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a 45 = A 41 + A 43 


(3.4.70) 


A 46 = A 42 + A 44 


(3.4.71) 


Additional relations required to complete the set of equations are given by 
Equation of State 


P = PT 


Turbulent Viscosity Relations 


e 1 = 1 + T e/p 

e 2 = 1 + r*e/y 


- i + r 


Pr 

?v n 


e 

v 


(3.4.72) 


(3.4.73) 


(3.4.74) 


(3.4.75) 


where T is the intermit tency factor and <t> is the ratio of eddy viscosity coeffi- 
cients if the turbulence is not isotropic. The calculations of these terms 
depends on the turbulence model which is described in Section 3.8. The total 
velocity is given by 


2 2 2 
U T = Uj + U 2 +2 


12 

H 1 H 2 


U 1 U 2 


(3.4.76) 


and the temperature ratio is 


given by 


0 = 


H t - U T 2 /2 

H t - U t 2 /2 
i e x e 


(3.4.77) 
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3.5 Properties of Boundary Layer Equations 


The set of equations derived in Section 3.4 consists of four coupled 
partial differential equations; three of which are second order and one first 
order. The properties of these equations have been formally examined by Wang 
(Ref. 18). Wang's analysis shows that these equations are parabolic and there- 
fore may be solved by a forward marching algorithm. However it is noted that 
there exist two sets of characterst ics (see Fig. 4); a set of characteristics 
which consists of surfaces normal to the wall, and a set of surfaces consisting 
of the stream surfaces. The set of surfaces normal to the wall are the three 
dimensional equivalent to the characteristic lines normal to the wall found in 
the two dimensional boundary layer equations. Thus one may expect that wall and 
edge boundary conditions must be applied similar to those applied to the two 
dimensional boundary layer equations. The second set of characteristc surfaces, 
consisting of stream surfaces, produce the hyperbolic like properties of the 
three dimensional boundary layer equations. From these surfaces, zones of 
influence and zones of dependence can be determined as shown by Wang (Ref. 18). 
These zones are determined by constructing a stream surface from the envelope of 
streamlines passing through a line perpendicular to the surface at the point in 
question as shown on Fig. 4. The zone of dependence is given by that volume of 
space sweeped by all vertical lines passing through the stream surface upstream 
of the line in question. Likewise the zone of influence is given by that volume 
of space sweeped by all vertical lines passing through the stream surface down- 
stream of the line in question. Since information is propagated along the 
characteristic surfaces, the solution for the flow properties on the line in 
question depends on the flow properties upstream in the zone of dependence. 
Likewise the flow properties on the line in question influence the flow down- 
stream in the zone of influence. Clearly in a forward marching algorithm, the 
solution of the flow properties must be known in the zone of dependence in order 
to calculate the flow on any vertical line. Thus the initial conditions for the 
problem are all the flow properties on a surface not a stream surface. These 
initial conditions are called inflow conditions because the solution can be 
obtained only downstream along stream surfaces. Thus, as pointed out by Blottner 
(Ref. 16), a unique solution of the three dimensional boundary layer equations 
requires specification of the inflow conditions along any inflow surface and 
specification of boundary conditions similar to similar to those employed for the 
two dimensional boundary layer equations. The boundary conditions are treated in 
Section 3.6 and the inflow conditions in Section 3.7. 


3.6 Boundary Conditions 

The boundary conditions are the same as those for the two dimensional 
compressible boundary layer equations. The no slip condition at the wall takes 
the form 
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U 1 

(x p x 2 . 

0) = 0 

(3.6.1) 



U 2 

(Xp x 2 . 

0) = 0 

(3.6.2) 

With normal 

injection at 

the wall 

we have 





U 3 

(Xp x 2 , 

0) = u w 

(3.6.3) 

Where Uy is 

the surface 

inject ion 

velocity 

which has been set to zero 

in the 


present work. For the energy equation, one may specify either the wall 
temperature or the heat flux at the wall. Thus we have 


or 


T (Xj, X 2 , 0) = T w 


(3.6.4) 


3t 



(3.6.5) 


At the boundary layer edge, we have: 


lim Uj = U le 
*3 » 


(3.6.7) 


lim U, - U 
X 3 - - 


2e 


(3.6.8) 


lim H T - H Te 
*3 + » 


(3.6.9) 


In the transformed variables, these boundary conditions take the form 


F Up 5 2 * 0) = 0 


(3.6.10) 
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or 


G (? 1# Z 2 , 0) = 0 


H 1 H 2 


v (Cp Z 2 > °) /2 ?i t 77“ p w u 3w 


(3.6.11) 

(3.6.12) 


H (Cp z 2 > o) = 


y t 


w 


•y-l H 


(3.6.13) 


Te 


3H y Pr 

_ (5l , tjf 0 ) 


/2Ei 1 


^W 


(3.6.14) 


Y_1 p e U le ^ Re H 1 L W p e H Te Po^coCptref 


F = G = H = 1 lim 5 3 


(3.6.15) 


Although specification of the boundary conditions Eqs. (3.6.10) through 
(3.6.15) is sufficient to solve the boundary layer equations, they are not suffi- 
cient to solve the problem because the thermodynamic state of the flow is not 
uniquely specified. This requires that the density and pressure be specified at 
the edge of the boundary layer. Thus of a complete specification of the problem 
requires the specification of U^ e , l^g, P e , H Tg , p e at the edge of the boundary 
layer in addition to the above boundary conditions. These quantities are deduced 
in the present study from a solution of the surface Euler equations subject to 
the imposed experimental pressure distribution. This analysis is discussed in 
detail in Section 3.11. 


3.7 Inflow Conditions 


According to Section 3.5 , the inflow conditions as well as the boundary 
conditions given in Sect. 3.6 must be specified. These inflow conditions include 
Uj, U 2 , P, T, p, H Te along any inflow plane as shown on Fig. 4. These properties 
may be specified as input data or constructed from analytic relations. In 
general this data is not known. An alternative approach is to construct the 
inflow conditions from certain degenerate solutions of the partial differential 
equations themselves. These degenerate solutions are called similarity solutions 
and are of two types. The first type, called local similarity solutions, are 
obtained by reducing the partial differential equations in three independent 
variables to partial differential equations in two independent variables. Thus, 
as an example, by setting % 1 = 0 , and solving a set of partial differential 
equations with £ 2 and 5 3 as independent variables, the inflow conditions along 
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the = 0 plane is obtained. The second type of degenerate solutions is 
obtained by reducing the partial differential equations to ordinary differential 
equations. These similarity solutions may be called true similarity solutions 
with one subset the similarity solutions for the two dimensional boundary layer 
equations. A third special cases does exist for a boundary which is a plane of 
symmetry for the flow field. In this last case, the plane of symmetry is a 
characteristic surface and requires special treatment as will be shown below. 

Local Similarity Solutions 

For local similarity solutions, the inflow conditions are only a function 
of and along a 5 2 boundary and only a function of C 2 anc * £ 3 along a 5^ 
boundary. An inspection of the differential equations in Section 3.4 indicates 
that these solutions can be obtained by setting 


A^ = 0 along a ^ = const, boundary 


A 2 = 0 along a = const, boundary (3.7.1) 

and using the local edge conditions. 

True Similarity Solutions 

One subset of true similarity solutions is for two dimensional 
incompressible laminar flow. For this subset it is assumed that: 


= HjdXj 

(3.7.2) 

= h 2 x 2 

(3.7.3) 

12 = 0 

(3.7.4) 

= U leo S i m 

(3.7.5) 

2e = 0 

(3.7.6) 

: = H TeO 

(3.7.7) 
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Note that these solutions require orthogonality of coordinates at the inflow 
plane and require a power law expansion of the velocity edge conditions. With 
these assumptions the differential equations of Section 3.4 , reduce to a set of 
ordinary differential equations where all coefficients A n = 0 except 


2m 

Ao = = $ 

J m+1 


2m 

A| -j = = 8 


4 17 


*28 


m+1 

2m 

m+1 


= 8 


A 36 " 1 


A 44 “ 1 


A 46 = 1 


(3.7.8) 


The equations of Section 3.4 reduce to the well known Falkner-Skan equations of 
which two solutions are of special interest. 


8=0 Flat Plate Solution 
8=1 Stagnation Point Solution 


(3.7.9) 


Although these true similarity solutions were not used in the present study, they 
have general usefulness in a variety of real problems. As an example on a 
cascade of blades with no spanwise pressure gradient, the Falkner-Skan solution 
with 8 = 1 can start the calculation at the leading edge stagnation point. 

Plane of Symmetry Solution 

Since no flow crosses a plane of symmetry, a plane of symetry is a 
characteristic surface and the solution of the equations are indeterminate. 
Blottner (Ref. 19) has derived the special relations needed to resolve the 
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indeterminancy in the equations. Across a plane of symmetry (X 2 - 0) the 
following conditions apply: 

ft = 0 (3.7.10) 

P (Xp x 2 ) = P (Xp - x 2 ) 

T (Xp X 2 ) = T (Xp - X 2 ) 

(3.7.11) 

Uj (Xj x 2 ) = Uj(Xp - x 2 ) 
u 2 (Xp x 2 ) = -u 2 (Xp - x 2 ) 

Clearly from Eq. (3.7.11) 

U 2 (Xp 0) = 0 (3.7.12) 

In addition to the flow variables having symmetry, the coordinates, the metrics, 
and the surface should have symmetry. 


'12 


= 0 


(3.7.13) 


3H 1 

3X^ 



(3.7.14) 


Then following the method of Blottner (Ref. 19), the solution is expanded in a 
power series in X 2 and keeping only the first term results in. 


8U 


2e 


U 


ref 


2e 


(3.7.15) 


3X 


2 ' o 


where the subscript (o) refers to the first term of the Taylor series expansion, 
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such that 


lim 3 U 3 ) = 1 

+ 00 


(3.7.16) 


Then by inspection of the equations in Section 3.4, the nonzero coefficients in 
the boundary layer equations are given by 


A 3 = 


A 1 - 2 5l 

2 <4 1 9U le 

^ U le 9 ^1 


A 12 = A 3 


A 14 A 3 


l 16 


h 

9 

“ 

' 9U 2e ' 

j 

% 

9 x i ] 

9X 2 



A 17 “ A 16 


v 18 


- h H i 


H, 1 


H 2 U le 


9U 


2e 


9X,, 


A 19 = A 18 


A 27 = A 19 


A 28 " A 17 


A 30 " A 16 + A 18 


(3.7.17) 
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= 2Cj 1 3H Te 


A 41 A 18 

A 44 “ 1 


A 45 " A 41 

A 46 A 42 + A 44 


(3.7.18) 


Note that the expansion of removes the indeterminancy in the coefficients 

A 16 and A 18 * 


3.8 Turbulence Model 

The turbulence model used by Vatsa (Ref. 1) is based on the Cebeci and Smith 
model (Ref. 20) developed for two dimensional boundary layers. This model may be 
extended to three dimensional boundary layer flow by replacing the streamwise 
velocity with the total velocity. Thus in the inner region we have 



(3.8.1) 


where the derivative of the total 


velocity and the mixing length are given by 


3u.p 

9x^ 



2Gi2 8uj 

h l h 2 9x 3 



(3.8.2) 
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£ 


m 


KX3 


1 - exp 



(3.8.3) 


where k is the Von Karman constant generally taken to be 0.4. The coefficient in 
the damping factor is given by 


where 



(3.8.4) 


* 


u 


(t w / Pw ) 


1/2 


(3.8.5) 



(3.8.6) 


In the outer region, the eddy viscosity is given by the Clauser model in the form 


e 

y 



(3.8.7) 


where £ Q is a representative length scale and y is the Clauser constant generally 
taken to be 0.016. The outer layer length scale, as developed by Mellor and 
Herring (Ref. 21), is taken to be 


00 


£ 


o 


/ 


x 3 


3u 

8x3 


dx- 


(3.8.8) 


The intermittency factor T is used to model laminar /turbulent transitional 
flow and varies between the limits 0.0 and 1.0. Since a transition model does 
not exist for three dimensional boundary layers, the simplest forced transition 
model for two dimensional flow is used. This model is that developed by Dhawan 
and Narasimha (Ref. 22) which is given by 
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(3.8.9) 


where S T ^ and are the specified locations for the beginning and end of 
transition. Finally we note that three dimensional turbulent boundary layers may 
have nonisotropic turbulence. Thus the eddy viscosity in the cross flow 
direction is multiplied by a factor 41 which may take a value as low as 0.40. 


3.9 Finite Difference Equations 

The governing equations described in Section 3.4 are a set of four coupled 
nonlinear partial differential equations. These equations are reduced by Vatsa 
(Ref. 2) to a set of four coupled nonlinear finite difference equations using the 
finite difference operators given below. 


8Q K K 

JjT~ = (Q I,J " Q I-1 


(3.9.1) 


3Q K K 

= (Q I,J * Q I,J-1 )/(5 2,J " ^.J-l* U 2 > 0 0.9.2) 


(Q I-1,J+1 " Q I-1 ,J^^2 ,J+1 ^2 , .P U 2 < 0 (3.9.3) 


3Q K+l K-l K+l K-l 

— = (Q I}J - Q i> j)/(C3 - £3 ) 


(3.9.4) 


H- 


3Q 


9 ?3 J 


e K ♦ e K+1 
£* +1 _ 5 K-l 


0 K+1 - 0 K 

r K+l c K 
$3 ~ £3 


e K ♦ e K_1 


rK+l _ r K-l 

^3 *3 


0 K - 0 K_1 
2LxJ ^i.J 

e K . F K-l 
?3 ^3 


(3.9.5) 
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where Q represents a general solution vector (F, G, V, H) and the i, j, k sub- 
scripts refer to the streamwise ^ , spanwise £ 2 , and normal C 3 locations 
respectively. Note that upwind difference is used on the crossflow convective 
derivatives depending on the sign of the crossflow velocity so as to honor the 
region of dependence criteria. When the crossflow velocity is negative, the 
differencing is explicit and there is a stability criteria on the step size which 
must be satisified. 


< _ H i**i 
U 2 h 2 ax 2 


(3.9.6) 


The resulting nonlinear difference equations are quasi-linearized to form a set 
of linear matrix equations which are solved in an iterative fashion (see Vatsa, 
Refs. 1 and 2). The matrix equations are solved by a block substitution algori- 
thm described in detail by Vatsa (Ref. 2). 


3.10 Generalized Surface Coordinates 

The computational coordinates are surface coordinates that wrap around the 
blade surface such that when X 3 = 0.0 , any point on the surface is a function of 
only X 1 and X 2 . Furthermore the computational domain shall be arbitrarily 
bounded by 


0.0 £ X 1 £ 1.0 
0.0 < X 2 < 1.0 


(3.10.1) 


without any loss of generality and as assumed in Section 3.1, the coordinate X 3 
shall be normal to the surface and hence normal to both X 1 and X 2 . The coordi- 
nates of the physical surface ( Yj, Y 2 , Y 3 ) are generally known only in Cartesian 
coordinates. Therefore we shall use the Cartesian coordinates as the basis for 
the transformation. (Note subscripts and superscripts are used to denote covar- 
iant and contravariant tensor components respectively.) The purpose of this 
section is to find a general tensor transformation from the computational coordi- 

nates (X) to the physical (Cartesian) coordinates (Y) for an arbitrary surface. 

It shall be noted at the outset that the Cartesian coordinates of the surface are 
generally not known in analytical form but only in numerical form so that a 
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general surface spline fitting algorithm is required to construct the transforma- 
tion. Thus we are searching for a transformation of the form 


Y = Y (X) (3.10.2) 


having the above properties and such that the Jacobian 


J 



* 0 


(3.10.3) 


within the computational domain. A nonorthogonal coordinate system having these 
properties can always be constructed provided that the surface projects uniquely 
to a plane surface. Thus if one projects all points on the surface onto the Y 3 = 
0.0 plane such that Eq. (3.10.3) is satisfied, then uniqueness is assured. Under 

■f 

these conditions, (see Fig. 2), the Cartesian coordinates of the surface (Y) can 
be written in parametric form using only the computational coordinates X^, X^). 
The transformation, Eq . (3.10.2), then takes the special form 


Y x = Y l (X 1 , X 2 ) 

(3.10.4) 

Y 2 = Y 2 (X 1 , X 2 ) 

(3.10.5) 

Y 3 = Y 3 [y^x 1 , X 2 ), Y 2 (X 1 , X 2 )] 

(3.10.6) 


It is noted that if the transformations Eq. (3.10.4) and (3.10.5) are known, 
then the equation of the surface in the physical coordinates is known in terms 
of the computational coordinates by direct substitution via Eq. (3.10.6). The 
task then is to determine the transformations given by Eqs . (3.10.4) and 
(3.10.5). From Fig. 2, it is noted that the four boundaries projected onto the 
Y 3 = 0.0 plane correspond to 


X^ = 0.0 side 1 
X 2 = 0.0 side 4 
X 1 = 1.0 side 3 
X 2 = 1.0 side 2 


(3.10.7) 
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Thus, as an example, the Cartesian coordinates of the side 1 boundary can be 
written in terms of only X 2 as a parameter. Since these boundaries are known 
only at descrete points in the form of numerical data, the side 1 boundary can 
parametized with X 2 as a variable using the following procedure for N boundary 
points . 

Y 1 <»■ x2 > ' Y 1 I x2 (I >l S ide 1 ‘ Y l (I >»id e 1 

y 2 x2 > " y 2 [ x2(I) l S ide 1 ‘ Y 2 (I) S ide 2 U.10.8) 

X 2 (I) = ( I— 1 ) / ( N— 1 ) I - 1, N 

Thus to each pair of points Y^I), Y 2 (I) there corresponds a unique X 2 (I). A 
similar proceedure is followed on the remaining three sides. With the four 
boundaries written in parametric form with the appropriate X 1 or X 2 as a 
variable, any interior point is uniquely calculated using the transfinite mapping 
of Gordon and Thiel (Ref. 5). 

Y i (X 1 X 2 ) = (1-X 1 ) Y i (0, X 2 ) + X 1 Y £ (1, X 2 ) 

+ (1 - X 2 ) Y i (X 1 0) + X Z Y i (X 1 1) 

- (1-X 1 ) (1-X 2 ) Y i (0, 0)- (l-X 1 )X 2 Y i (0,l) (3.10.9) 

- X 1 (l-X 2 ) Y l (1,0) - X X X 2 Y^l.l) 


where i = 1,2. Then given any X 1 and X 2 , the Cartesian coordinates (Yj, Y 2 , Y 3 ) 
are determined uniquely using Eqs . (3.10.9) and Eq. (3.10.6). The covariant 
metric tensor components are then determined from 



3Y k _^k 
3X i 3Xj 


(3.10.10) 
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Hi = /G n 


(3.10.11) 


H2 “ / G22 


H 3 = 1 


|g| = det ] Gi j j (3.10.12) 

where the contravarient tensor notation is dropped for Eqs. (3.10.10) through 
(3.10.12). We note the the covariant metric tensor components (see Eq. 

(3.10.10)) require first derivatives and also that the equations of motion 
requires derivatives of of the metric tensor (see Section 3.2). Therefore first 
and second derivatives of the the transformation Eq. (3.10.2) are required. 
Interpolation and differentiation of the boundary curves is done using a four 
point Lagrangian interpolation formula. Interpolation and differentiation of Eq. 
(3.10.6) is done using a surface spline fit developed by McCartin (Ref. 23) and 
Spath (Ref. 24). 

In addition to the metric tensor components G^j, the boundary layer equa- 
tions also require components of the coordinate rotation vector ft (see Section 
3.2). In addition it is usefull to be able to resolve any vector in the 
computational coordinates to the corresponding vector in the Cartesian 
coordinates. Therefore the direction cosines from the computational to Cartesian 
coordinates must be determined. In the present work these direction cosines are 
used to resolve Cartesian components of the boundary layer edge velocities and 
the wall shear force vector. A vector in Cartesian space is given by 

dR - e^Y 1 + e 2 dY 2 + e 3 dY 3 (3.10.13) 

where is the unit Cartesian vector. However on the surface, from Eq. 

(3.10.6) 


dY 3 = 


3Y- 


3Y J 


dY 1 + 


3Y- 


dY" 


3Y^ 


(3.10.14) 


hence Eq. (3.10.13) becomes 
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dR = ejdY 1 + e 2 dY z + e 3 


On the surface X 3 = 0.0 and 


9Y- 


9Y J 


dY 1 + 


9Y- 


dY 


9Y^ 


(3.10.15) 


dY 1 = 


dY 2 = 


9Y 


9X 


9Y^ 


9X- 1 


dX 1 + 


dX 


1 + 


9Y 


9X 


9Y 


9X" 


dX 


dX^ 


(3.10.16) 


Substituting Eq . (3.10.16) into Eq . (3.10.15), we have 


->■ 

dR = 


+ 

e, 


9Y 


9X 


9Y 


-► 

+ 


■> 

+ e 0 


9Y" 


9X 


9Y 


-»• 

+ e. 


+ 

+ e -i 


9X 


9X^ 


j 9Y 3 
\9Y X 

/ 9Y 3 
\ 9Y 1 


9Y 1 

9Y 3 

9Y 2 \ 

+ 

| 

9X 1 

2 

9Y Z 

9 X 1 J 

9Y 1 

3 

9Y J 

9Y 2 \ 



9X 2 

9Y 2 

9X 2 /. 


!dX x 


dX 2 (3.10.17) 


The contravariant components of dR can also be written as 

dR = a^dX 1 + a 2 dX z + a 3 (0) 


(3.10.18) 


*»■ O 

with a as the basis vector and noting that dX = 0 on the surface. Hence by 
comparing terms between Eqs . (3.10.17) and (3.10.18) we have 
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where the distinction between covariant and contravariant vectors may be dropped, 
The direction cosines are then given by 


Y ij 



-► 



(3.10.21) 


Then if is any vector in the X^ direction in the computational coordinates and 
Vj is any vector in the Cartesian coordinates, then 
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(3.10.22) 


3.11 Surface Euler Equations 

The solution of the boundary layer equations requires specification of the 
edge conditions (U le , U 2e , H Tg , p g , P e ) as shown in Section 3.5. These edge 
conditions can be obtained directly from experimental data or they can be 
obtained from solutions of the Euler equations. An alternative method is to 
solve the surface Euler equations using a known surface static pressure 
distribution. This method was outlined by Cebeci et al. (Ref. 7) for application 
to aircraft wings. Gleyzes and Coustiex (Ref. 8) developed a similar method for 
incompressible flow over fusiform bodies. This section extends these methods to 
obtain a solution of the compressible surface Euler equations for application to 
arbitrary three dimensional surface such as turbine blades. 

The surface Euler equations can be obtained from the three dimensional boun- 
dary layer equations given in Section 3.2 by taking the limit as •+• <, °. Thus at 
the edge of the boundary layer normal derivatives vanish and hence the momentum 
and energy equations are written as follows: 


Momentum 


U le 3 Ut „ 8U 


h, ax 


le . 2e le M n 2 + it it + a -.it 
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h 2 ax 2 
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(3.11.1) 


X 2 Momentum 
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(3.11.2) 
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Energy 
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The coefficients in these equations are given by 
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The static pressure has been expressed in terms of the pressure coefficient 
defined as 


C 


P 



2 (P - P ) 
00 


(3.11.4) 


since experimental data is normally given in this form. 

Inflow Conditions 

The solution of these equations requires specification of the velocity com- 
ponents U| e and U 2 e as well as the thermodynamic variables P g , T e , p e along the 
inflow boundaries. In general, only the local surface static pressure coeffi- 
cient and the freestream or reference stagnation conditions P Too and T Too are 
known. If the local flow angle can be estimated, then the local edge conditions 
can be calculated by expanding the flow isentropically from the known stagnation 
conditions to the local measured static pressure. 

A useful relation for this purpose is one between the pressure coefficient 
and velocity. For incompressible flow this is given by the Bernoulli equation. 
Since on a rotor the total pressure and total temperature are not constant for 
the whole flow field, they are assumed constant along a surface of constant 
radius. (The Bernoulli surfaces are cylindrical surfaces). Thus P^ and T T are a 
function of radius and we have 


P + - U e • U e (3.11.5) 

2 
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where 


•* •+• 

u e * u e = ^e ' U g + 2 U 0 • ( Q x R) + (f2R) 2 (3.11.6) 

If the vorticity is aligned with the Y 1 axis and the radius vector is aligned 
with the Y 2 axis, then 

U e • (fix R) = (fiR) ^Y 13 U le + Y 2 3 u 2ej (3.11.7) 

Likewise the reference conditions are given by 


P T <» = p « 


■> 

U_ 


(3.11.8) 


U • U = U +2 (Rfi) U. + (C1R) 2 

oo 00 oo v w OO (4)00 v / 


(3.11.9) 


Combining these equations results in an explicit relation between and U £ . 



a o U le U 2e + 


2e 


d l U le 


d 2 U 2e 


1 ” C p + d 0 + d 4 + <^5 (3.11. 10) 


where it is noted that p = p„ - 1 and U w = 1 for incompressible flow. For com- 
pressible flow, the same reasoning applies. Thus total pressure and total tem- 
perature are constant along cylindrical surfaces and 


Y 
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The energy equation can be written as 



(3.11.12) 


and combining the equations we have, 


u le + U 2e + 2a o U 2e + d l U le + d 2 U 2e " d 8 + d 9 


(3.11.13) 


where the coefficients are given by 
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d 9 = (J1R) 2 
d o “ 2 <P T - P T J 


a o = 2G 12 /(H 1 H 2 ) 


Auxiliary Relations 

Three additional relations are required to complete the set of equations. These 
are the equation of state 


P = p 

K e e 


(3.11.14) 


and the equations for total enthalpy and rothalpy 


Y 

H ^- — ^*~r 


(3.11.15) 


h = H Te 


(3.11.16) 


Solution Algorithm 

These equations are hyperbolic equations in which the characteristics are 
the streamlines. Therefore these equations have the same properties as the boun- 
dary layer equations discussed in Section 3.5 and are solved by specifying the 
inflow conditions on the inflow boundary. These inflow conditions are given by 
(Ui e , U 2e , H^ e , P g ). Since pressure (i.e., Cp) is a known input, the velocity 
components Uj e and U 2e may be determined by specifying either U 2g or the flow 
angle on the inflow plane and using Eq . (3.11.10) or (3.11.13) to determine Uj g . 
With the velocities known, the temperature is determined from Eq.(3.11.14) and 
Eq . (3 . 11 . 15) . Finally the density is determined from Eq . ( 3 . 1 1 . 14) . 

These equations have the same properties as the boundary layer equations and 
the finite difference equations are formed in the same manner using the Eqs . 
(3.9.1), (3.9.2) and (3.9.3). This results in a set of coupled algebraic equa- 
tions . 
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b l U le 2 + b 2 U le + b 3 U le U 2e + b 4 U 2e + b 5 U 2 e 2 + b 6 = 0 (3.11.17) 

b 7 U le 2 + b 8 U le + b 9 U le U 2e = b 10 U 2e + b ll U 2 e 2 + b 12 = 0 (3.11.18) 

e 4 I e - e 5 (3.11.19) 

where the coefficients are given by 
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where the coefficients on the left side of the page are for U 2e >0.0 and those 
on the right are for U 2 g ^ 0.0. Eqs. (3 . 1 1 . 17 ) » (3 . 11 . 18) , (3 . 1 1 . 19) are solved 
iteratively in a nested loops by successive substitution. The outer loop solves 
for the density. If the flow is incompressible, the density p = 1.0. The inner 
iteration loop solves for the velocity components and rothalpy. Assume U 2e is 
known, then Eq. (3.11.17) can be solved for U le . With U^ e known, Eq. (3.11.18) 
can be solved for U 2e> Since both equations are quadradics, they are solved in 
the following manner; 
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and 


♦j - (b 2 * b 3 U 2e ) 2 - 4b | (b 6 * b A 0 2e ♦ bjll 2e 2 ) (3.11.20) 


% - < b l 0 + b 9 D le )2 - 4 b ll (b 12 * b 8 U le * b 7 “l. > (3.11.21) 
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(3.11.23) 


where the positive roots are appropriate. When the velocity components are 
known, the rothalpy is determined from Eq. (3.11.3), and the remaining roperties 
from Eqs . ( 3 . 1 1 . 14 ) , ( 3 . 1 1 . 1 5 ) , and (3.11.16) so that the density is updated in 
the outer iteration loop. 
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4.0 RESULTS AND DISCUSSION 


4.1 Introductory Discussion 

In this assessment of the applicability of three dimensional boundary layer 
theory to predict the flow streamlines and heat transfer in a gas turbine passage 
four regions of the flow field are examined. These regions are the turbine 
pressure surface, the turbine endwall surface, the turbine suction surface, and 
the rotor pressure surface as indicated by the shaded areas on Fig. 1. The first 
three cases test the analysis in a stationary coordinate system and the fourth in 
a rotating coordinate system. Experimental data for pressure, endwall, and 
suction surfaces were obtained by Graziani et.al. (Ref. 11) in a large gas 
turbine cascade which simulated a turbine rotor. Detailed wall static pressure 
distributions, three component velocity traverses, wall heat transfer, and wall 
limiting stream lines were presented. Data for the rotating blade case was 
obtained by Dring And Joslyn (Refs. 12 and 13). This data includes wall static 
pressure distributions and wall limiting streamlines. In addition, radial 
traverses of the total pressure and flow angle in the stationary reference frame 
upstream of the rotor were given. This data was sufficient to calculate the 
cases presented in this section. 

Initial results using this three dimensional boundary layer analysis have 
been presented by Vatsa (Refs. 1 and 2) wherein the boundary layer edge condi- 
tions were obtained directly from the velocity traverses. In the assessment 
presented in this section the boundary layer edge conditions were obtained by 
integrating the surface Euler equations using the measured static pressure 
distributions. The overall analytical procedure, which was the same for all 
cases, is shown on the flow chart shown on Fig. 3. Spline smoothed Cartesian 
coordinates of the turbine blade surface are used to calculate a coordinate 
system using the geometry analysis presented in Section 3.10. With the 
coordinates known, the experimental pressure distribution was used to calculate 
the boundary layer edge conditions using the surface Euler analysis given in 
Section 3.11. Finally, the three dimensional boundary layer equations, given in 
Section 3.4, were solved. The inflow conditions were estimated using the local 
similarity approximation, described in Secton 3.7, along all inflow planes that 
are not characteristic surfaces. For the cascade pressure surface case, the 
plane of symmetry is a characteristic surface and the inflow conditions were 
estimated using the plane of symmetry analysis also described in Section 3.7. 
Integral properties of the boundary layer such as displacement thickness and 
momentum thickness (defined for a three dimensional boundary layer) were 
calculated for all cases. These integral properties of the boundary layer were 
obtained from the velocity vectors resolved in an intrinsic coordinate system 
defined by the free stream flow direction and the cross flow direction (see 
Appendix). In addition the two components of the wall friction coefficient 
(streamwise and crossflow), heat transfer (Stanton number), and wall limiting 
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streamline skew angles were calculated. The Stanton number was calculated in two 
steps. In the first step, the adiabatic wall temperature was calculated by 
assuming zero heat flux at the wall. In the second step, a constant heat flux 
equal to that applied experimentally was used, and the wall temperature 
calculated. Thus the calculation procedure used for determining the Stanton 
number was the same as that used experimentally in Ref. 11. 


4.2 Turbine Cascade Pressure Surface 

The turbine cascade pressure surface developes a nearly two dimensional 
boundary layer since the passage vortex is near the suction surface as shown 
on Fig. 1. Since the spanwise static pressure distributiuon is nearly consant, 
little spanwise flow or skewing is expected. For this reason, the primary purpose 
of this test case is to demonstrate a calculation started along a plane of symme- 
try boundary. The computational domain extends from a point just down stream 
from the leading edge of the turbine blade to a point just upstream of the trail- 
ing edge of the blade and spanwise from the plane of symmetry to the endwall. 

The coordinates calculated for this case consists of a 40 by 40 grid with 100 
grid points normal to the surface grid of which every other coordinate line is 
plotted on Fig. 5. The calculated boundary layer edge conditions are shown on 
Fig. 6 in the form of velocity vectors projected on to the = 0.0 plane. These 
velocity vectors show, as expected, very little spanwise flow at the edge of the 
boundary layer. Since the pressure distribution was not accurately known near 
the leading edge, the boundary layer calculation was started with a finite bound- 
ary layer thickness using the local similarity approximation by specifying a 
finite distance (£j > 0) for the virtual origin (see Sect. 3.7). The boundary 
layer was assumed to be completely turbulent throughout. A plot of the wall 
limiting streamlines is shown in Fig. 7 in the form of the wall stress vectors 
projected on to the Y 3 = 0.0 plane. Again this vector plot shows very little 
spanwise flow. In addition, the skew angle (the angle between the boundary layer 
edge flow direction and wall streamline direction) is not more than ten degrees. 
These results compare favorably with the experimental results for flow direction 
and flow skewing presented by Graziani et . al . in Ref. 11. 


4.3 Turbine Cascade Endwall Surface 

The endwall surface boundary layer has a large cross flow produced by the 
strong static pressure gradient from the pressure surface to the suction surface 
as shown in Fig. 1. In addition this boundary layer has a complex pattern 
produced by a separation and attachment line and a saddle point. Downstream of 
this region the boundary layer is rapidly accelerated and the boundary layer 
growth is more systematic. In this downstream region it may be possible to use 
three dimensional boundary layer theory. The computational domain and skewed 
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coordinate system used for this case is shown on Fig. 8. It consists of a 40 by 
40 grid with 100 grid points normal to the surface grid which extends from a 
point just downstream from the saddle point to just upstream of the trailing edge 
to the suction surface of the blade passage. The boundary layer edge conditions 
were calculated using the wall static pressure distribution for the thick bound- 
ary layer case presented by Graziani et. al. (Ref. 11). Along the inflow bound- 
aries for the surface Euler calculation, it was assumed that the boundary layer 
edge flow was tangent to the coordinate line so that U 2 e was assumed zero. 
Integration of the surface Euler equations showed that futher downstream small 
negative cross flow velocities developed indicating that the boundary layer edge 
flow curvature was less than the coordinate curvature. This result was consis- 
tent with the velocity traverse data used by Vatsa (Refs. 1 and 2). Since this 
small negative U 2 e velocity component produced stability problems (see Section 
3.9) which could not be resolved even with a much finer grid, the calculation for 
the edge velocity conditions was run assuming that U 2 e was zero in agreement with 
the procedure of Vatsa (Ref. 2). A plot of the edge velocity vectors, which are 
tangent to the coordinate lines, is shown on Fig. 9. 

The boundary layer calculation was started with inflow conditions obtained 
from local similarity solutions of the boundary layer equations. A comparison of 
the calculated wall shear vectors with the measured surface streamlines is shown 
on Fig. 10. As seen in this figure, the flow inclinations are in qualitative 
agreement with the measured data except in the neighborhood of the saddle point. 
It is also noted that the calculation presented here using the measured wall 
static pressures to derive the boundary layer edge conditions shows much better 
agreement for the predicted flow angle than the calculations of Vatsa (Ref. 2) 
which used measured velocities for the edge conditions. This result indicates 
the difficulty in determining the edge of the boundary layer in a three dimen- 
sional flow field. Improved results perhaps may be obtained with better defini- 
tion of the flow conditions along the upstream inflow plane. A comparison of the 
measured and calculated heat transfer (Stanton number) distributions is shown in 
Fig. 11 in the form of contour plots. Again qualitative agreement is obtained 
except in the region of the saddle point. The case presented by Vatsa (Ref. 1) 
is the same as that presented here except that the edge velocities were deter- 
mined directly from three component velocity traverses whereas in the present 
case the edge velocities were calculated from the wall static pressure distribu- 
tion using the surface Euler equations. Of interest then is a comparison of the 
predictions of heat transfer using the two methods. This is shown on Fig. 12 in 
the form of crossection plots of Stanton number vs. spanwise distance a several 
axial stations. Good qualitative agreement is obtained between the two methods 
where the difference between the results illustrates the difficulty in obtaining 
precise boundary layer edge conditions from measured data in a complex three 
dimensional flow field. 
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4.4 Turbine Cascade Suction Surface 

The suction surface boundary layer of a gas turbine blade develops strong 
three dimensional effects due to the presence of the passage vortex which sweeps 
the flow up from the end wall surface. Since the boundary starts at the stagna- 
tion point near the leading edge of the blade, the leading edge inflow conditions 
are well established. However the flow from the end wall inflow surface produces 
a separation line which divides the midspan flow from the end wall flow. This 
inflow is less well established. The computational domain and coordinate system 
for this case is shown on Fig. 13. It consists of a 40 by 40 grid with 100 grid 
points normal to the surface grid extending from just downstream from the leading 
edge of the blade to just upstream of the trailing edge of the blade and from the 
end wall to the mid plane or plane of symmetry. The boundary layer edge condi- 
tions were calculated assuming zero cross flow velocities at the leading edge 
inflow plane. Along the endwall inflow plane, the flow angle and hence cross flow 
^2e velocit y was estimated from the measured surface streamlines given by 
Graziani et al. (Ref. 11). A plot of the calculated edge velocity vectors is 
shown on Fig. 14. Since the cascade surface is a developable surface, the 
vector plot is shown using arc length as coordinates. This vector plot clearly 
shows the effect of the passage vortex on the flow at the boundary layer edge. 
Since the surface Euler equations are hyperbolic equations with streamlines as 
the characteristic lines, this case illustrates difficulties which may be 
encountered as the flow approaches the plane of symmetry. The overall flow 
direction is from the endwall to the plane of symmetry. Therefore the flow along 
the plane of symmetry is a result of the calculation not an input boundary or 
inflow condition. Along this plane of symmetry the cross flow velocity is zero. 
Hence the static pressure must be such as to reduce the cross flow velocity to 
zero at the plane of symmetry. In general without sufficient data this may be 
difficult to achieve as was found in the present case near the blade trailing 
edge. 


The three dimensional boundary layer equations were solved using local 
similarity inflow conditions as described in Section 3.7 and the local calculated 
edge conditions. Since this is a transition boundary layer, the begining and end 
of transition must be estimated because the present analysis lacks a transition 
model for three dimensional boundary layers. The begining and end of transition 
were obtained by calculating the quasi-two dimensional flow along the plane of 
symmetry using the intermittancy function of Dhawan and Narashima (Ref. 22) and 
comparing the results with the experimental data for heat transfer. These 
transition points were then held constant across the span of the blade. A com- 
parison of the calculated wall shear vectors with the measured wall limiting 
streamlines is shown on Fig. 15. The flow direction at the wall is seen to be 
calculated quite well and the separation line separation line dividing the lead- 
ing edge flow from the end wall flow is also predicted quite well. A comparison 
of the measured and calculated Stanton number distributions is shown on Fig. 16. 


54 


R85-956834 


Again encouraging agreement is obtained both as to magnitude and overall pattern. 
The calculation, however appears to underpredict the peak heat transfer which 
occurs near the endwall inflow plane. Crossection plots of the Stanton number 
versus spanwise distance at different axial stations are shown on Fig. 17. These 
plots show that the average level of heat transfer across the span is predicted 
reasonably well, but that the high and low regions are underpredicted. Calcu- 
lated and measured heat transfer along the midplane of the blade is shown on Fig. 
18. This figure also shows predictions made from a two dimensional calculation 
using the present three dimensional boundary layer analysis and also predictions 
made by Sharma (Ref. 25). These comparisons show that although the flow is quasi 
two dimensional along the plane of symmetry (i.e. = 0.0), the heat transfer 

is not predicted accurately. Sharma (Ref. 25), with a quasi 3-D calculation, 
showed that this difference can be accounted for by the streamline convergence of 
the flow toward the plane of symmetry. Thus it may be expected that the differ- 
ence between the calculated 3-D heat transfer and the measured heat transfer can 
in large part be accounted for by errors in the prediction of streamline conver- 
gence. Vatsa (Ref. 1) obtained better predictions of the heat transfer along the 
midplane using measured velocities for the edge conditions. Since no pressure 
distribution data was available from the quarter-plane to the midplane for this 
case, the streamline convergence of the flow towards the midplane could not be 
estimated by the surface Euler analysis as accurately as otherwise might be 
possible. 

The three dimensional boundary layer analysis also calculates the wall fric- 
tion coefficients and integral parameters such as displacement and momentum 
thickness. Since a nonorthogonal coordinate system may be rather arbitrary and 
not reflect the flow properties, these parameters are calculated in a locally 
intrinsic coordinate system (see Appendix). This locally intrinsic coordinate 
system is referenced to the boundary layer edge flow direction and an orthogonal 
cross flow flow direction. All velocities and shear stresses are locally resolved 
to this coordinate system. Thus two displacement thicknesses are defined. One 
represents the velocity defect in the streamwise direction and the other the 
velocity defect in the cross flow direction. Likewise the wall friction coeffi- 
cients represents the wall shear in the streamwise and cross flow direction. 
Contour plots of the streamwise component of wall friction coefficient and 
displacement thickness are shown on Figs. 19 and 20. 


4.5 Turbine Rotor Pressure Surface 

The flow over the rotor pressure surface is strongly influenced by the rela- 
tive vorticity generated by rotation of coordinates. Thus the boundary layer 
flow is driven by two applied forces; the pressure gradient, and the Coriolis 
force. Because of rotation rotation, the pressure surface sees a strong radial 
pressure gradient which tends to drive the flow radially inward. The Coriolis 
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force may tend to drive the flow radially inward or outward depending on the sign 
of the blade rotation vector normal to the surface fig. The balance of these 
forces tends to drive the flow radially outward near the leading edge and axially 
near the trailing edge. The computational domain and coordinate grid for this 
problem is shown on Fig. 21. It consists of a 40 by 40 grid with 100 grid points 
normal to the surface grid extending from just downstream from the leading edge 
stagnation line to just upstream of the trailing edge and from the rotor hub to 
the rotor tip. The boundary layer edge velocity vectors calculated from the sur- 
face Euler equations are shown on Fig. 22. It was found that no solution of the 
surface Euler equations exist along the leading edge inflow boundary for the 
crossflow velocity U2 e - 0.0. Therefore on both the leading edge and hub inflow 
boundaries, the crossflow velocity was estimated from the surface streamlines 
given by Dring and Joslyn (Refs. 12 and 13). The edge velocities show the 
behavior of the suction surface flow observed by Dring and Joslyn. 

The boundary layer equations were started with local similarity inflow 
conditions. A comparison of the calculated wall stress vectors with the wall 
streamlines obtained using the ammonia-ozalid process is shown on Fig. 23. The 
calculated vectors clearly show the general overall flow pattern quite well. 
Calculated flow angles are compared with the flow angles measured from the ozalid 
paper on Fig. 24. Two curves are shown. The first curve (solid line) is the 
calculated flow angle at each point along the calculated streamline. The second 
curve (dotted line) is the calculated flow angle along the midplane. A compari- 
son of the edge velocity vectors with the wall vectors indicates that the flow 
skewing is up to about ten degrees. These figures clearly show that the flow 
streamlines are predicted fairly accurately and that two dimensional boundary 
layer theory applied in a strip fashion would not be adequate to calculate the 
heat transfer on the turbine blade. Calculated Stanton number distributions for 
a fully turbulent boundary layer are shown on Fig. 25. This shows higher heat 
transfer rates near the hub and trailing edge of the rotor blade. 


56 



R85-956834 


5.0 CONCLUDING REMARKS 


An assessment has been made of the applicability of a three dimensional 
boundary layer analysis for the calculation of heat transfer, total pressure 
loss, and streamline flow patterns on the surfaces of both stationary and 
rotating turbine passages. On the pressure and suction surface of a turbine cas- 
cade, the 3-D boundary layer analysis along with the coordinate generator and 
surface Euler analysis produce predictions of heat transfer and surface stream- 
line patterns which are in reasonably good agreement with available experimental 
data. Experimental data for wall friction coefficient was not availble to 
evaluate the prediction of loss; however, because of the close relation between 
wall friction and heat transfer one may also expect good predictions of loss. 
Assessment of the 3-D boundary layer analysis for the predictions of the flow on 
a turbine endwall is complicated by the uncertainty in the inflow conditions 
downstream from the separation and attachment lines. However the results 
obtained in this study are very encouraging as they indicate that the 3-D bound- 
ary layer analysis is applicable to the endwall surface as well. On the pressure 
surface of a rotating turbine blade, prediction of the flow streamlines and flow 
skewing agree reasonably well with the available data. Evaluation of the analy- 
sis for the prediction of heat transfer and loss on a rotating turbine blade must 
await the collection of experimental data. 

The results obtained in this investigation reinforce earlier conclusions 
that two dimensional boundary layer theory applied in a strip fashion can not 
predict heat transfer and loss because of the strong three dimensional flow 
effects on a turbine blade. On the suction surface of a blade, these three 
dimensional effects are produced by the passage vortex and on the pressure 
surface of a rotating blade these effects are produced by the radial pressure 
gradient and Coriolis force. In view of the current state of Navier-Stokes analy- 
sis, three dimensional boundary layer theory could provide a useful tool for 
local enhancement of of the flow field for the prediction of heat transfer, loss, 
and flow streamline pattern on turbine blade surfaces. The results of this study 
strongly indicate that the local details of the flow field on the turbine blade 
surface can be accurately calculated using a two step procedure. First the over- 
all flow field can be calculated on a coarse grid using a Navier-Stokes or 
parabolized Navier-Stokes analysis. Second the three dimensional boundary layer 
analysis can be used for local enhancement to provide details of the flow field 
on a much finer grid. 
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8.0 LIST OF SYMBOLS 


a„ , c_ 
n ’ n 


A , K 
n ’ n 


c p» c v 


g ij> G ij 

G 

h, H 
h , H 

f Hip 

H 

i , I 

k 

i 

L 

M 

-► 

n i 

P , P 


Coefficients of boundary edge equations 

Coefficients of boundary layer equations 

Specific heats 

Wall friction coefficient 

Pressure coefficient 

Unit vectors Cartesian coordinates 


X 2 velocity ratio 

(u l /u le ) 

Covarient metric tensor components 

<8 i3 /* 2 > 

X velocity ratio 

(u 2 /u ref ) 

Metric scale coefficients 

(h./£) 

Static enthalpy 

(h/u£) 

Total enthalpy 

(h T /u£) 

Total enthalpy ratio 

( h'p/h'j’g ) 

Rothalpy 

(i/u£) 

Thermal conductivity 


Reference length 


Levy - Lees parameter 

(py/p e W e ) 

Mach number 


Unit vectors surface coordinates 


Static pressure 

0 

(p/PooUoo) 
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Px> 

Total pressure 

2 

( Pj/ Poo u oo 

Pr, Pr t 

Prandtl number laminar/turbulent 

(Cp p/k) 

q 

Levy - Lees length scale parameter 


• 

q 

Heat flux 


r , R 

Radius 

(r/A) 

Re 

Reynolds number 

(p^U^l/P. 

<% 

Gas constant 


s i’ S i 

Arc length distance 

(s./4) 

S t 

Stanton number 


t , T 

Static temperature 

( t/t re f ) 

u i’ U i 

Velocity components surface coor. 

(uj/u.) 

u z> U z 

Axial velocity 

(u z /uj 

u 

S3 

u 

0 

Radial velocity 

(u r / Uco) 

u <t» u <f> 

Tangential velocity 

(u^/Uoo) 

u rfef U ref 

X~ reference velocity u- or u 0 „ 

2 J le 2e 

(u ref/u<» 

v i> V i 

Velocity components Cartesian coor. 

(v i /u 00 ) 

V 

Xg reduced velocity 


V B’ V B 

Rotor velocity 

(Vg/u*,) 

X i> X i 

Surface coordinates 

(x^/A) 

yi* Y i 

Cartesian coordinates 

( yi /A) 

h 

Relative inlet flow angle 


y 

Ratio of specific heats 
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J-J 

-1 

ij 


e ref 

0 

0 

U 

^ref 

«i 

P, P 


b), 


Direction cosine tensor components 
Inverse tensor components 
Displacement thickness 
Eddy viscosity 

Viscosity ratio turbulent/laminar 
Reference viscosity ratio 
Momentum thickness 

Static temperature ratio (t/t f) 

Molecular viscosity 

Reference molecular viscosity v(t re f) 

Levy - Lees transformed coordinates 
Density (p/p^) 

Time 

Stress components 

Rotor speed (iiit/uj 


Subscripts 

c Cross flow direction 

e Edge conditions 

h Hub 

L rotor inlet conditions 

S Streamwise direction 

t Tip 

T Total or stagnation conditions 

« Free stream or reference, conditions 
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Superscripts 

~ Absolute rotor reference frame 

Tensor Notation 

i,j,k Subscripts covarient tensor components 

i,j,k Superscripts contravarient tensor components 
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Fig. 1 Flow in a Turbine Passage 
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Fig. 3 Flow Chart of Geometry, Surface Euler and 3-D Boundary Layer Analysis 
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Fig. 4 Zone of Influence and Zone of Dependance for Three Dimensional Boundary Layers 
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Fig. 5 Coordinate System Calculated for Cascade Pressure Surface 
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Fig. 7 Boundary Layer Wall Shear Vectors Projected on Y 3 Plane Calculated with the 
3-Dimensional Boundary Layer Analysis for the Cascade Pressure Surface 
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Fig. 9 Boundary Layer Edge Velocity Vectors Calculated from Experimental Surface Pressure 
Distributions for Cascade Endwall Surface 
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Fig. 10 Comparison of Measured and Calculated Endwall Limiting Streamlines 
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Fig. 12 Comparison of Stanton Number Distributions on Endwall Surface 
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Fig. 13 Coordinate System Calculated for Cascade Suction Surface 
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Fig. 14 Boundary Layer Edge Velocity Vectors Calculated from Experimental 
Surface Pressure Distribution for Cascade Suction Surface 
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Fig. 15 Measured and Calculated Wall Limiting Streamline on Cascade 
Suction Surface 
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Fig. 18 Calculated and Measured Stanton Number on Mid Plane of Cascade Suction Surface 
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Fig. 19 Streamwise Component of Wall Friction Coefficient on Cascade 
Suction Surface 
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Fig. 21 Coordinate System Calculated for Turbine Rotor Pressure Surface 
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Fig. 24 Calculated and Measured Flow Angle on Pressure Surface of a Rotating Turbine Blade 
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Fig. 25 Calculated Stanton Number Distribution on Pressure Surface of a Rotating 
Turbine Blade 
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10.0 APPENDIX - BOUNDARY LAYER PARAMETERS 


The boundary layer parameters such as momentum thickness, displacement 

thickness, and wall shear are calculated in an intrinsic coordinate system (see 

Ref. 26). This intrinsic coordinate system is oriented with respect to the 

boundary layer edge flow direction 0 e and the crossflow direction defined as the 

tfo X tL direction. These directions are labeled the streamwise and crossflow 
j e 

direction respectively. The resolution of vectors in this coordinate system is 
given in Ref. 26 and the following boundary layer parameters are defined. 

Displacement Thickness 



6*/£ = / [ 1 - ] dYo 

S Pe U e 

(A. 1) 


6*/£ = / [ ] dY 3 

C Pe U e 

(A. 2) 

Momentum Thickness 

pU U 

9.Jl = [ 1 " 1 dY, 

88 Pe U e U e 

(A. 3) 


6 Jl = -/ ^ ^ dY 3 
“ ’ pu c », 3 

(A. 4) 

Wall Friction Coefficient 




C £s ‘ 2 

(A. 5) 


C £c " 2 W<P e U e 2 ) 

(A. 6) 

Stanton Number 

s 

c p» c p (t a - t w ) 

(A. 7) 
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